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Abstract 

In this work we obtain the exact analytical scattering solutions of a particle (electron or hole) in a semiconductor 
double heterojunction — potential well / barrier — where the effective mass of the particle varies with position 
inside the heterojunctions. It is observed that the spatial dependence of mass within the well / barrier introduces 
a nonlinear component in the plane wave solutions of the continuum states. Additionally, the transmission 
coefficient is found to increase with increasing energy, finally approaching unity, whereas the reflection coefficient 
follows the reverse trend and goes to zero. 
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I. INTRODUCTION 

Recent developments in the nanofabrication of 
semiconductor devices have given a thrust to the 
study of quantum mechanical systems with posi- 
tion dependent effective mass. The spatial depen- 
dence on the effective mass of the particle arises 
due to its interaction with an ensemble of particles 
within the device, as the particle propagates from 
left to right. This so-called position-dependent 
effective mass (pdem) formalism becomes an es- 
sential ingredient in describing the electronic and 
transport properties of quantum wells and quan- 
tum dots, impurities in crystals, He-clusters, quan- 
tum liquids, semiconductor heterostructures, etc. 
P-Q. When two such materials (having pdem) 
with different bandgaps are placed adjacent to each 
other to form a heterojunction, the effective mass 
approximation is valid within each material. If, 
for example, a thin layer of a narrower (wider) 
bandgap material is sandwiched between two lay- 
ers of a wider (narrower) bandgap material, they 
form a double heterojunction. If the intermedi- 
ate layer is sufficiently thin for quantum proper- 
ties to be exhibited, then the alignment is called a 
single quantum well (barrier). A typical quantum 
well structure may be composed of a semiconduc- 
tor thin film, embedded between two semi-infinite 
semiconductor materials, say GaAs / Al x Ga\- x As, 
where x denotes the mole fraction. As the mole 
fraction varies along the z-axis, so does the effec- 
tive mass of the charge carrier (electron or hole). 
To have a complete understanding of such a quan- 
tum system, quantum mechanics requires knowl- 
edge of both bound and scattering states. 

Various attempts have been made over the years 
to study the scattering states in such position- 
dependent effective mass systems [9hfl6l|. In two 
of the recent works [13, EH, the authors stud- 



ied the effect of hard wall confinement and lateral 
dimensions on low temperature transport proper- 
ties of long diffuse channels in InSb/ Irii^ x Al x Sb 
heterostructures rjj, and resonant inelastic x- 
ray scattering in LaAlO^/ SrTiO^ heterostruc- 
tures [la] . In another work, ballistic carrier emis- 
sion was studied with GaAs / Al x Ga\^ x As as a 
model system [19( . In one of the relatively earlier 
works, the authors considered the simple model 
of a square potential well, with the effective mass 
varying with position inside the well [l3j |. However, 
as the additional term introduced by the changing 
mass is small compared with the original potential 
Vo and does not change the shape of the potential 
significantly, this term was neglected while finding 
the solutions. 

In the present work our aim is to obtain the 
exact analytical solutions for the scattering states 
of a particle inside a single quantum well/barrier. 
It may be mentioned that exact analytical solu- 
tions play an important role in conceptual un- 
derstanding of physics. They provide a valuable 
platform for checking and improving approximate 
models and numerical results. Herein lies the mo- 
tivation for obtaining exact analytical solutions of 
wave equations with pdem, especially because of 
the wide range of applications of these solutions 
in various areas of material science and condensed 
matter [20j. In the double heterojunction consid- 
ered here, we assume the intermediate layer to be 
a potential well / barrier of the form 



V(z) , ai < z < a.2 
P = I Voi = V(ai) , -co < z < ai 
Vq2 = V(a 2 ) , a 2 < z < oo 



(1) 



where ai and a 2 represent the heterojunctions. 
The mass of the charge carrier is assumed to be 
spatially varying inside the potential well / barrier 
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a\ < z < et2, but constant outside, viz, 

m(z) , eti < z < a,2 

mi = m(ai) , — oo < z < a\ (2) 
TO2 = m[a,2) , a-2 < z < oo 

Thus both the potential V(z) and mass m(z) are 
real functions of the configuration space coordi- 
nates, and are taken to be continuous throughout 
the semiconductor device. The work done in refs. 
[IH [20( | deserve special mention here. In ref. [l5| 
approximate analytical solutions were derived for 
any arbitrary potential and arbitrary mass func- 
tion, whereas in ref. [20| . some special forms of 
the mass function were considered for oscillator, 
Coulomb and Morse potentials to produce exact 
analytical results. However, our study differs sig- 
nificantly from both these works — the mass func- 
tions considered in our case, as well as the ap- 
proach used, are different from refs. [l5| and (2G |; 
the present article not only gives exact analytic re- 
sults but also plots the transmission and reflection 
coefficients. More importantly, unlike [IEH^j this 
work deals with a double heterojunction. 

The article is organized as follows : For the 
sake of completeness, the position-dependent-mass 
Schrodinger equation is introduced in Section II, 
and the method of obtaining the solutions is dis- 
cussed. A couple of explicit models are studied in 
Section III. To give a better insight into the phys- 
ical nature of the problem, the potential and mass 
functions are plotted as a function of z (in Fig. 1 
and Fig. 3) along with the scattering solutions (in 
Fig. 2 and Fig. 4). The transmission and reflec- 
tion coefficients are also calculated and the same 
are plotted as a function of the energy of the par- 
ticle, in Fig. 5 and Fig. 6 respectively. Section IV 
is kept for Conclusions and Discussions. 



II. THEORY 

We start with the basic one dimensional time in- 
dependent Schrodinger equation associated with a 
particle endowed with pdem in the intermediate 
region within the heterojunctions : 

H EM (z)iP(z) = [T EM (z) + V(z)} 1>(z) = E^{z) 

where the kinetic energy term Tem is given by 

aa 



Tem — ^ (m a prn^pm 1 + m 7 pm^pm") 



1 f 1 



(4) 



with p — —ih— being the momentum operator, 
dz 

and the ambiguity parameters a , /3 , 7 obeying 
the von Roos constraint |j| 



a + /3 + 7 



-1 



(5) 



There is neither a unique nor a universal choice for 
the ambiguity parameters, and several suggestions 
exist in literature — e.g. a = 7 = , /? = — 1 [2~i| . 
a = 7 = -1/2 , /? = 0[H, a = 7 = -1/4 , /3 = 
-1/2 [H, = 7 = -1/2 ,a = O0, etc. How- 
ever, for continuity conditions at the abrupt inter- 
faces, one should consider a = 7, or else one gets 
the unphysical result of the wave function vanish- 
ing at the heterojunctions: additionally the ground 
state energy also diverges [H,[2^|. We shall restrict 
ourselves to the BenDaniel-Duke choice for the am- 
biguity parameters, viz., a = 7 = , /3 = — 1. In- 
cidentally, this particular choice consistently pro- 
duces the best fit to experimental results [l?} ■ Fur- 
thermore, we shall work in units h = c = 1, and 
use prime to denote differentiation w.r.t. z. Thus, 
inside the potential well/barrier a\ < z < 02, the 
Hamiltonian for the particle with pdem reduces to 

Hi 



H = 



1 d 2 ( 1 V d 



2m(z) dz 2 \2m(z) J dz 



+ V(z) (6) 



whereas, outside the well/barrier, z < a\ and z > 
a2, the particle obeys the conventional Schrodinger 
equation : 



+ ^01,02 > ip(z) = E^{z) (7) 



1 d 2 
2mi,2 dz 2 



having plane wave solutions. In case we consider 
a wave incident from left, the solutions in the two 
regions are 



^ L {z) = e iklZ + Re- lklZ , -00 < z < o x 



(8) 



0,2 < Z < OO 



where R and T denote the reflection and transmis- 
sion amplitudes, and 



h,2 = \ 2mi,2 (E - V ifi 2 ) 



(9) 



To find the solution in the region a\ < z < a?, we 
make use of the following transformations [291 ] 

^ m - {2m(z)} 1/4 , p= { ^2m(z)dz (10) 



which reduce the Schrodinger equation for 
position-dependent mass, to one for constant mass, 



+ {np)-^=o 



(ii) 
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with 

~, s , s 7 to' 2 m" , s 
V ^ = V ^ + 32^-S^ (12) 

We are interested in studying the scattering states 
of a particle in a double heterojunction, formed by 
dissimilar materials, where the mass of the parti- 
cle varies with position inside the well / barrier. 
Such heterojunctions can be described by a mate- 
rial potential which derives from the difference in 
bandgaps @ . Crystal potential of multiple hetero- 
junction can also be described in this manner. For 
this purpose, we look for some definite practical 
forms of V(z) and m(z) which will give exact an- 
alytical solutions of (jTTJ) . We illustrate this with 
the help of a couple of explicit models in the next 
section. 



III. EXPLICIT MODELS 

A. Case 1 : Potential well with position de- 
pendent mass 



so that after some straightforward algebra V(p) in 
eq. ([12]) reduces to 

?M - ^ -(>- w) xA % (16) 

Thus equation (1111) can be written as 

+ |^ 2 + A(A-l)sech 2 ^|0 = O (17) 

where k 2 = E \-r (18) 

4/3 2 

and the parameter A depends on the constants p 
and /3, through the equation 

A(A-1)= M 2 -^ (19) 

For existence of bound states A > 1; hence, | p | > 

— . This gives the permissible values of A as 
zp 

A 4 + W i+4/i2 -^ (20) 



The one dimensional finite square well is one of the 
simplest confinement potentials. In the first ex- 
ample, we consider the region between the abrupt 
heterojunctions (— ao < z < do) to be a symmetric 
potential well (more realistic than the square well) 
with the following ansatz (p being some constant) 



V(z) = I 



1 + z 2 

,2 



-a < z < a 



l + al 



= V , 



-oo < z < — ao , ao < z < oo 



(13) 

This particular model resembles the profile of a 
diffused quantum well, with the advantage of exact 
analytical solutions. In case the real situation is 
slightly different from this model, one can apply 
approximation methods like perturbation theory, 
etc., to obtain the solutions. 
Let the mass of the particle be 



m(z) 



P 2 



2 (1 + z 2 ) 

P 2 



-ao < z < ao 



mo , 



2(1 + a 2 ) 

— oo < z < — ao , ao < z < oo 



(14) 

where /? is some constant parameter. For the spa- 
tial mass dependence given by eq. (|T4)) . eq. (fTO]) 
transforms the coordinate z to 



p = (3 sinh z 



(15) 



Black (dolled) line -» m (z) 
Red (dashed ) line -* V Izl 



FIG. 1: Colour online : Plot showing m(z) and V(z) 
w.r.t. z 

For a better understanding of the mass dependence 
and the potential in the semiconductor device, we 
plot m(z) and V(z) as a function of z in Fig. 1, for 
a suitable set of parameter values, viz., f3 = 4, p = 
3, a = 2. 

Let us introduce a new variable 

2 P 



y = cosh 



P 



(21) 



and write the solutions of (Tr7|) as 

4> = V k u(y) (22) 

In terms of the new variable y, equation (|17[) re- 
duces to the hypergeometric equation 

, d 2 u f / 1 \ , , , I du 

v(i-»)^ + {(A + 5 J-(A+i)y}^ 



dy 2 

i{A 2 + K 2 /3 2 }w = 



(23) 



3 



with complete solution [30( 

« = p 2F 1 ( K aA\■^-v}+Q{l-y) ^/2 ■ 
' i , 13 
2 F X ( a +-,6+~,-;l-^ 

(24) 

where P and Q are constants, and the parameters 
a and 6 are as defined below : 

a=^(X + iK(3) , b=-(\-iKP) (25) 

Thus for the even solution Q — 0, whereas for the 
odd solution P = 0. After some straightforward al- 
gebra, the final solution to the position-dependent 
mass Schrodinger equation within the poten- 
tial well — eto < z < ao, is obtained as 



d 2 



2(1 



1/4 



(l + z 2 ) 



A/2 



P 2 F 1 la, 6, -; -z' 

( 113 
+ iQz 2 Fi ( a + -,&+ -,-;-z 2 

(26) 

whereas outside the well (z < — ao , z > ao), the 
solutions are given by eq (|3|), with k\ — k 2 . 
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FIG. 2: Colour online : A plot of Re ip(z) vs z; Dashed 
(red) lines show the abrupt heterojunctions 



For scattering states, k 2 should be positive, imply- 
ing E > — — Now, because of the spatial depen- 

dence of the mass function, the boundary condi- 
tions need to be modified — 

the functions ip( z ) an d 



should be con- 



m(z) dz 

tinuous at each heterojunction ±ao |21j,|31|. These 
conditions enable us to evaluate the reflection and 
transmission amplitudes R and T respectively. 
The complete solutions for the scattering states in 
the entire region — oo < z < oo are plotted in Fig. 
2, for the same set of parameter values as in Fig. 



1, viz., /3 = 4,/i = 3,a = 2,P = 40. The 
solutions show a definite nonlinear character inside 
the well (—ao < z < ao), where the particle mass 
to is a function of its position z. Thus the effect 
of the position-dependent mass potential well (in 
this particular model) is to introduce a non linear 
component in the otherwise plane wave solutions. 



B. Case 2 : Potential barrier with position 
dependent mass 

As the second example, we consider the region 
within the double heterojunction (a\ < z < a 2 ) 
to be represented by an inverted Morse potential 
(barrier) (H-HI : 



V(z) = 



V e az (2 - e az ) , o x < z < a 2 

V e aai (2 - e QQl ) = V 01 , -oo < z < <ii 
Vo e aa2 (2 - e QQ2 ) = V 02 , a 2 < z < oo 

with positive Vq- The transmission coefficient of a 
potential barrier has wide applications in nuclear 
fission, heavy- ion fusion, tunnelling in solids [Hj], 
etc; hence we shall calculate both the transmis- 
sion as well as reflection coefficients and plot them 
as a function of the energy of the particle. The 
Morse barrier potential is particularly useful in in- 
vestigating the anharmonicities of the vibrational 
spectra in molecular and nuclear physics 33] . For 
the mass function of the particle, we consider 



m(z) 



rrioa 2 e~ 2az , a\ < z < a 2 

moct 2 e~ 2aai = 777oi ) — oo < z < a% 
moa 2 e~ 2aa2 = 77702 , a 2 < z < 00 



(28) 



FIG. 3: Colour online : Plot showing m(z) and V(z) 
w.r.t. z , The dashed lines show the abrupt hetero- 
junctions 



In Fig. 3, we show the plot of m(z) and V(z) for 
this particular model, for a suitable set of param- 
eter values, viz., toq = 0.4, Vo = 5, E = 33, with 
the heterojunctions at a± = —0.8 , a 2 = 0.8. 
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For the spatial mass dependence given by eq. 
cq. (|T0|) transforms the coordinate z to 



(29) 



so that after some straightforward algebra V(p) in 
cq. (fT2|) reduces to the simple form 



V 



2Vbv / 2^ ) 2m Vb - 3/4 



(30) 



Thus the Schrodinger equation for constant mass 
([TT]) takes the final form 



2 2moVb-3/4 2^oV2W\ , n 

, i i K H 3 1 <P = 



7' 



dp 2 IP" P 

(31) 

with k 2 — E. To solve equation (1511 given above, 
let us introduce a new variable 



y = —2inp = in2\/2moe az 



(32) 



in terms of which equation (|31[) gets simplified to 
the form of a Whittaker differential equation [35| 



d 2 ^ 
dy 2 



1 . Ag + 1/4 
4 



zA 2 



= o 



with 

\\ = 2m V - 1 



The solutions of eq (|33|) are given as [3£ 
= e ± ^ 2 y ±iAl M(a ± ,& ± ;y) 



(33) 



(34) 



(35) 



where M (a^ ,b ;y) are the Whittaker functions 
and 



1 



± i\± — 1X2 



l±2«Ai (36) 



Thus for the entire semiconductor device, the com- 
plete solution for the scattering states in the dif- 
ferent regions are given by 



ik x z 



Re tklZ , -00 < z < ai 

(2m ) 1/4 VSe^ [p ie iy^M (a+,b + ;y) 
+ P 2 e-%y- lXl M (a', b~\ J/)} , a x < z < 

tp R = Te lk2Z , a 2 < z < 00 

(37) 

where the constants Pi , P2 and the reflection and 
transmission amplitudes R and T respectively, are 
determined by matching the boundary conditions 
(for pdem systems) at the heterojunctions, and y 
and p are as defined above. 



The complete solutions for the entire region — 00 < 
z < 00 are plotted in Fig. 4, for a suitable set of 
parameter values, viz., mo = 0.4, Vq = 5, E = 
33, a\ — —1.5, a,2 — 1.5. The solutions show 
a definite nonlinear character inside the barrier 
(ai < z < a 2 ), where the particle mass is de- 
pendent on its position. Thus the effect of the 
position-dependent mass barrier (similar to that of 
the position dependent potential well in the previ- 
ous example) is to introduce a non linear compo- 
nent in the plane wave solutions. 




FIG. 4: Colour online : A plot of Re ip(z) vs z; Dashed 
(red) lines show the abrupt heterojunctions 



IV. CONCLUSIONS AND DISCUSSIONS 

To conclude, we obtained the exact analytical solu- 
tions for the scattering states of a particle (electron 
or hole) inside a semiconductor device with a dou- 
ble heterojunction, when the mass of the particle is 
assumed to be dependent on its position inside the 
heterojunctions, but constant outside. We studied 
two explicit models in this work — one a pdem 
diffused potential well, the other a pdem potential 
barrier (Morse barrier). In each case it is observed 
that the effect of the spatial dependence on the 
particle mass is to introduce a non linear compo- 
nent in the otherwise plane wave solutions. 
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FIG. 5: Plot of |7f vs E 
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FIG. 6: Plot of |7?| 2 vs E 



We also calculated the transmission and reflec- 
tion coefficients, |T| 2 and |i?| 2 respectively, for the 
potentials studied here. These are plotted in Fig. 5 
and Fig. 6 respectively, as a function of the energy 
E for the pdem particle in a Morse barrier. As the 
energy increases, the transmission coefficient also 
increases, finally reaching unity, whereas the reflec- 
tion coefficient follows the reverse trend and goes 
to zero. This observation is similar to that in ref. 
[36| , where the authors show that the transmission 
coefficient for the one-dimensional scattering prob- 
lem with pdem normally tends to unity as energy 
goes to infinity, provided the mass is a continuous 
function of position. 

This simple, yet straightforward approach is 
just a way of understanding basic physics of the 
electronic properties of a semiconductor device, 



comprising of a double heterojunction, where the 
intermediate layer is sufficiently thin for quantum 
properties to be exhibited. It is expected that the 
observations made in this work will provide some 
useful insight in studies related to electron trans- 
port in semiconductor heterostructures, i.e. in the 
physical properties of such materials. Actual ma- 
terials are made up of a large number of atomic 
potentials. Nevertheless, the crystal potential may 
be approximated by a single potential — the global 
average of the individual potentials, and this ap- 
proach would still be valid. However, for extremely 
thin intermediate layer, the individual potentials 
may become significant enough for this approxima- 
tion to break down. This calls for a more rigorous 
approach, and we propose to take up its study in 
the near future. 
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